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CN ! Abstract 
> 

' The prediction of backgrounds to new physics signals in topologies with large missing transverse 

cN ; 

■ energy and jets is important to new physics searches at the LHC. Following a CMS study, we 

. investigate theoretical issues in using measurements of 7+ 2-jet production to predict the irreducible 

o 

background to searches for missing energy plus two jets that originates from Z + 2-jet production 
where the Z boson decays to neutrinos. We compute ratios of 7 + 2-jet to Z + 2-jet production 



^ . cross sections and kinematic distributions at next-to- leading order in a^, as well as using a parton 

shower matched to leading-order matrix elements. We find that the ratios obtained in the two 
approximations are quite similar, making 7+ 2-jet production a theoretically reliable estimator for 
the missing energy plus two jets background. We employ a Frixione-style photon isolation, but we 
also show that for isolated prompt photon production at high transverse momentum the difference 
between this criterion and the standard cone isolation used by CMS is small. 

PACS numbers: 12.38.Bx, 13.85. Qk, 13.87.Ce 
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I. INTRODUCTION 



The LHC era is now upon us, and the hunt for the mechanism of electroweak symmetry 
breaking and new physics beyond the standard model is underway. Typical signatures for 
supersymmetry and many other new physics models include topologies with large missing 
transverse energy (MET) accompanying jets (MET J). The same signatures can easily be 
mimicked by Standard Model processes, such as the production of an electroweak boson 
decaying into neutrinos, in association with jets. For the discovery of new physics in the 
early running of the LHC, with only a few inverse femtobarns of integrated luminosity, it is 
important to understand the Standard Model backgrounds to METJ searches. 

Events containing a Z boson and jets, with the Z decaying into a neutrino pair (METZJ), 
constitute an irreducible background to the METJ signal. One can envisage using com- 
plementary approaches to understanding this and other backgrounds: a direct theoretical 
prediction; or data-driven approaches, which estimate the rate from measurements of other 
processes (or possibly from other kinematic regions in the same process). Data-driven tech- 
niques offer a powerful means of avoiding theoretical uncertainties in background predictions, 
as well as cancelling experimental systematics common to different processes. However, such 
methods can require theoretical assistance, in order to estimate the "translation" parameters 
from one process to another, and their inherent uncertainties. The question theorists can 
address is how stable is the ratio of two processes to various theoretical approximations. 

In the case of the process Z[^ z/z/) + jets, the most obvious choice of other process would 
be Z{-^ ^"'"£~) + jets, i.e. the production of a Z boson in association with jets, where the Z 
decays into a charged lepton pair. The production kinematics and dynamics of these two 
processes are identical, so no theoretical input about QCD is required, only knowledge of the 
Z boson branching ratios. Leptonically decaying Z bosons have the drawback, however, of 
offering very low event rates, less than a sixth of METZJ (per lepton channel), even before 
imposing lepton rapidity cuts. The paucity of statistics has led experimenters to examine 
using other processes to estimate METZJ rates and distributions. The CMS collaboration 
has studied l2| the use of W or photon production in association with jets for estimating 
the METZJ background. The production of a in association with jets offers an order of 
magnitude higher statistics than the leptonic Z process; the production of a prompt photon 
in association with jets, sixteenfold higher statistics than leptonic Z decays. In addition. 
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FIG. 1: Squark pair production illustrates a new physics process with the signature of two jets plus 
MET. Here each squarks decays to a quark and the lightest neutralino; the escaping neutralinos 
generate the missing transverse energy. 

W production suffers from contamination from ti events, and could well be affected by the 
same kinds of new physics the experiments are seeking in the METZJ sample. Reducing ti 
contamination requires selection cuts that enhance the photon channel's advantage. 

Use of either the W or j plus jets processes requires knowledge of their short- distance 
strong-interaction dynamics, as they probe different combinations of the parton distributions 
and somewhat different scales than Z plus jets. In the photon case, its masslessness also 
affects distributions in important ways and requires theoretical input for any comparison 
with massive boson production. The study of both processes is valuable, of course; the 
use of different processes allows for cross checks and also presumably different sensitivity to 
whatever new physics may be lurking in the data. 

In this paper, we study the 7+ 2-jet and 1^1^)+ 2-jet production processes. The latter 
process is a background to new physics, such as the squark pair-production process illustrated 
in fig. [H Our aim is to provide the necessary theoretical results, to next-to-leading order 
(NLO) accuracy in the strong coupling a^, for using 7+ 2-jet production to estimate Z+ 2-jet 
production. We use the same software tools as in our previous studies oiW + n- and Z + n- 



jet production 



-6], the BlackHat library I?, 8| along with AMEGIC-F-K [ol within the 
SHERPA lOj framework, to perform the leading-order (LO) and NLO calculations. We also 
present results for a parton-shower calculation matched to fixed-order LO matrix elements 

n □ 

(ME+PS) also within the SHERPA [10[ framework. A key issue is the theoretical 
uncertainty in the conversion from 7 to Z. We use the difference between the ME-I-PS 
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results and the NLO predictions to estimate the uncertainties for ratios. (The common 
variation of factorization and renormahzation scales in the numerator and denominator of 
these ratios produces quite small shifts in the ratios, which are likely to underestimate 
the uncertainties substantially.) The results presented here are being used by the CMS 
collaboration in their study of missing energy in association with three jets [21]. While our 
present study is for missing energy in association with two jets, we do not expect much 
difference in the theoretical uncertainties between the two- and three-jet cases. 

Photon isolation is essential for rejecting copious hadronic backgrounds. The type of 
photon isolation criterion affects the theoretical description of the photon production pro- 
cess. In the past, various types of isolation cones have been used, which limit the amount of 
hadronic energy near the photon candidate. Fixed isolation cones generally limit the total 
amount of energy in a cone, while the one proposed by Frixione 12|] consists of a set of 
energy constraints that become increasingly restrictive the closer one gets to the photon. 
The Frixione cone is theoretically attractive because it eliminates contributions from long- 
distance coUinear fragmentation of partons into photons. Although there is a perturbative 
factorization available for other types of cones, the required photon fragmentation func- 
tions [isl (non-perturbative functions analogous to the parton distribution functions) are 
not known particularly precisely. From an experimental point of view, some hadronic energy 
must be allowed everywhere within the cone in order to cope with the underlying event and 
with event pile-up. In our study, we adopt a modified cone criterion of the Frixione type. 
The CMS collaboration has recently published [l^ a measurement of the photon spectrum. 
We compare the two types of cone isolation to the data, and show that with our choice of 
parameters the difference between them is small for the kinematic region studied by CMS. 

In the following section we describe the details of our calculation and discuss the photon 
isolation criterion. Section IIIII presents our cross-check using isolated prompt-photon pro- 
duction. Section HV] discusses the cuts we use. In section |V] we present the ratios of Z + 2-jet 
to 7 + 2-jet rates for a variety of distributions. Our conclusions and outlook follow. 



II. THE CALCULATION 



In this section we 
setup used in refs. [3H5 



discuss our calculational setup. At NLO, we follow the same basic 
while for the ME+PS study we use the setup of ref. HI]. 
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(a) 





(b) 



FIG. 2: Sample virtual diagrams needed for (a) pp — )• Z{-^ uv) + 2-jet production and for (b) 
— )• 7 + 2-jet production. 

A. Matrix Elements and Integration 



We compute the cross sections at NLO using the Catani-Seymour dipole subtraction 
method 1_5|. This method requires the combination of several contributions: the LO term; 
virtual corrections from the interference of tree-level and one-loop amplitudes; the real- 
emission corrections with dipole subtraction terms; and the singular phase-space integrals 
of the dipole terms. 

We evaluate the required one- loop amplitudes using the BlackHat program library 
For the processes we are studying, we need the one-loop corrections to the following partonic 
processes, 



QQgg Z{-^ vv) or 7 , 
qqq'q' Z{-^ uu) or 7 . 



(2.1) 



where two of the four partons should be crossed into the final state, and the Z decay 
to neutrinos is folded in. Some sample diagrams for these processes are shown in fig. [2l 
illustrating the similarity of the Z and 7 cases. 



For the Z processes 



plitudes from refs. 
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the BlackHat code library [7,|8| implements analytic one- loop am- 



17 1 . (See also ref. |l8|.) The photonic amplitudes are implemented 



using analytic formulae in a similar way; they can be obtained as appropriate sums of color- 
ordered primitive amplitudes for the two-quark three-gluon 1^ and four-quark one-gluon 



processes 



2d ], as explained in ref. [21| and in the fourth appendix of ref. [19|. This conver- 
sion is possible because the trace-based color decomposition does not distinguish between 
different generators of U{N). Setting one of the generators equal to the identity matrix, and 
collecting the coefficients of the identical remaining color structures, generates the photonic 
amplitudes. This procedure removes the unwanted three-boson couplings of the photon 
present in each primitive amphtude. We omit the process gg -> ggj as it only contributes 
to 7 + 2-jet production at next-to-next-to-leading order; for the kinematics of interest here 
the gluon luminosity is not large enough for this process to be important. 

The NLO result also requires real-emission corrections to the LO process, which arise from 
tree-level amplitudes with one additional parton. We use the program AMEGIC-|--|- jslto 
compute these contributions, along with the Catani-Seymour dipole subtraction terms 
and their integrals over phase space. The SHERPA framework 10| includes AMEGICH — h 
and tools to analyze the results and construct a wide variety of distributions. We have 
previously validated the BlackHat-I-SHERPA framework for W + {n <2) jets against 
the MCFM code Q- 

We use SHERPA in a second mode, to provide a parton-shower prediction matched to 
tree-level matrix elements (ME+PS), also known as matrix-element-plus-truncated-shower. 
(Our parton-shower results do not include hadronization effects, but remain at the parton 
level.) The ME-I-PS event samples are produced following ref. using the COMIX 
matrix-element generator 23|]. This method combines LO hard matrix elements together 
with parton showers, which resum logarithmic corrections due to Bremsstrahlung effects. 
The parton shower employed to this end in SHERPA [2^ is based on Catani-Seymour 
dipole factorization 15|]. In contrast to earlier parton showers, the procedure inherently 
respects QCD soft color coherence. The procedure allows the unambiguous identification 
of a recoil partner for partons that are shifted off mass-shell in the splitting process (the 
"mother" partons). This eliminates one of the major sources of uncertainty in earlier schemes 
for parton evolution. As the observables presented below should be relatively insensitive to 
hadronization effects, ME+PS results are presented at the parton level. We match to matrix 
elements containing up to three final-state partons, and use 15 GeV for the merging cut. 
(Further details may be found in ref. |25|.) 
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We work to leading order in the electroweak coupling. The Z-boson couplings are as 
given in ref. [s]. In particular, the vv invariant mass is distributed in a relativistic Breit- 
Wigner of width Yz = 2.49 GeV about the Z boson mass of 91.1876 GeV. These values, 
along with those of aQEuiMz) = 1/128.802 and sin^ 6'^^ = 0.230 lead to a branching ratio 
for the neutrino mode in Z decay of Br(Z — t- uu) = 0.2007. We use MSTW2008 leading 
and next-to-leading order parton distribution functions, with the QCD coupling as chosen 
appropriately in each case. Our LHC results are for a center-of-mass energy of 7 TeV. 

The question of what scale should be chosen for the electromagnetic coupling in a prompt- 
photon calculation is a subtle one. Strictly speaking, this scale is undetermined in a leading- 
order computation; and although we are working to NLO in the QCD coupling, we are 
only working to LO in the electromagnetic one. From a practical point of view, however, 
it makes a difference whether we choose a scale of order the photon transverse momentum 
(and hence oem ~ «em(^z) ~ 1/128), or the zero-momentum-squared value (that is, 
aEM(O) = 1/137.036). 

Heuristically, we see that the sequence of fermion bubbles on the photon forms a gauge- 
invariant set. Moreover, because of QED Ward identities, this is the only set of diagrams 
which controls the coupling renormalization. In the idealized situation we are considering, 
the photon is always resolved, i.e. it is not allowed to split into low-mass lepton pairs. Then 
summing the bubbles leads to a zero-momentum-squared coupling for the emission of the 
hard photon. Other types of QED effects, such as emission of additional hard photons, 
change the kinematics and cannot be absorbed into a running coupling. This argument is 
confirmed by the analyses in refs. 26|, l27|] . We therefore take the electromagnetic coupling 
for the prompt-photon computation to be the zero-momentum-squared value, q;em(0) = 
1/137.036, and not the running value at typical collider energies. 



B. Infrared Safety and Photon Isolation 

An observable in perturbative QCD is infrared- and collinear-safe if it is unaltered under 
the emission of an arbitrarily soft gluon, or the splitting of a colored parton into a pair of 
colored partons (whether in the initial or final state). We follow ATLAS and CMS in using 
the infrared-safe anti-kx jet algorithm j28 |. 

For 7 -|- 2-jet production, there is an additional infrared issue beyond the infrared safety 
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of the jet algorithm. Because photons can arise from vr" decay and other hadronic sources, 
it is essential from an experimental point of view to insist that they be isolated from jets 
if we want to study photons that originate in short-distance physics ('prompt'). On the 
other hand, too strict a photon isolation criterion — such as requiring no hadrons or no 
hadronic energy in a cone around the photon, or using tracking alone to determine isolation 
— would be infrared- unsafe, and prevent the use of perturbative QCD as a theoretical tool. 
In order to navigate between these two competing requirements, experimental collaborations 
typically use a weighted isolation criterion, imposing a limit on the hadronic energy fraction 
in a cone around the photon, or simply on the total hadronic energy in the cone. 

In their recent measurement of the inclusive isolated prompt-photon spectrum 1^, the 
CMS collaboration required photon candidates to satisfy a set of requirements on nearby 
energy deposits measured via tracking, and in the electromagnetic and hadronic calorimeters. 
In their theoretical modeling using Pythia, they required photons to have less than 5 GeV of 
summed pt within an isolation cone of radius R = 0.4, where Ri^ = ^ {rji — rj^y + (0j — 0^)^. 
We will adopt this criterion as our 'reference' standard cone isolation in our discussion of 
the isolated prompt-photon spectrum. 

While such a criterion is infrared-safe with respect to the strong interactions — emis- 
sion of a soft gluon, or a colored parton splitting into two colored partons — it is not 
coUinear-safe with respect to QED: the cross-section receives contributions from coUinear 

radiation of photons off massless quarks. In a theoretical description, this singularity has 

I I 

to be factorized, and absorbed into parton-to-photon fragmentation functions [13], whose 
computation from first principles would require knowledge of non-perturbative physics. The 
factorization and the non-perturbative functions are the final-state analogs of the parton 
distribution functions. In practice, these functions are extracted from fits to experimental 
data, although these fits are not nearly as precise as those for the parton distributions, 
nor are sets surrounding a central fit available to estimate errors. (As is true for the par- 
ton distributions, the evolution of the fragmentation function with the hard scale can be 
determined in perturbative QCD.) 



However, the isolation criterion given above is not the only possible one. Frixione [12 1 
proposed a modified isolation requirement which suppresses the collinear region of the phase 
space and thereby eliminates the need for a fragmentation-function contribution. We follow 
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this proposal, requiring that each parton i within a distance Ri^ of the photon obey 

J2E,Te{6-R^y)<n{6), (2.2) 

i 

for all (5 < (5o, in a cone of fixed half-angle So around the photon axis. The restricting function 
"HiS) is chosen such that it vanishes as 5 — > and thus suppresses coUinear configurations, 
but allows soft radiation arbitrarily close to the photon. We adopt 



1 — cos 6 



n 



n{6) = E:^e[—-^\ , (2.3) 

y 1 — cos Oq J 

where Ej, is the photon transverse energy. 

An experimental analysis cannot adopt the Frixione prescription precisely, because of 
the finite resolution of detectors. One can imagine using a discretized version of it; this 
would still require a fragmentation contribution (corresponding to the innermost step of the 
discretized cone), but it would presumably have a much smaller one than the standard cone. 
One could imagine adjusting the parameters e, 6o and n so as to minimize that contribution. 
A preliminary study of the comparison between the Frixione prescription and a discretized 



version may be found in ref. 



29|. 



In our study, we will use the Frixione cone, with e = 0.025, 6o = 0.3 and n = 2. In 
practice, we find that our predictions are only weakly sensitive to these parameters. In the 
next section, we will compare our predictions using these parameters to predictions made 
using a standard cone for the isolated prompt-photon spectrum |l4| measured by CMS. As 
we shall see, the differences between the Frixione and standard-cone isolation prescriptions 
are not large, and are quite small in the large-pj. region that is of primary interest in the 
present study. We therefore conclude that it is reasonable to use the Frixione isolation to 
model the Z-to-photon ratio in association with two jets for CMS's analysis. 

In the definition of the cross section for 7 + m jets, at least m jets must lie outside the 
isolation cone, while we allow any number of additional jets m' to fall inside the cone of 
radius 6q = 0.3. In order to obtain the desired observable we apply the jet-finding algorithm 
to all partons except the photon in the event, to obtain m jets outside the isolation cone 
and m' jets inside the cone. Then we apply the cuts on transverse momenta and rapidity 
of the m jets. This procedure is infrared safe because the jet-finding algorithm is applied 
everywhere and not just outside the isolation cone. In practice, jets appear very rarely inside 
the photon isolation cone, because of the combination of the jet pt cut and the small energy 
fraction we use. We will hst the other cuts we apply in a later section. 



III. ISOLATED PROMPT-PHOTON PRODUCTION 

Recently, CMS has published a measurement [l^ of the inclusive isolated prompt-photon 
spectrum based on a data sample of 2.9 pb~^ from early running of the LHC They compared 
the measurement with NLO predictions from the publicly-available JetPhoX code 30|. 
These predictions make use of a standard photon isolation cone as described in the pre- 
vious section, and include fragmentation contributions. (The JetPhoX code implements 
a Frixione isolation cone as well as a standard isolation cone, although the former imple- 
mentation is not used in this study.) We can use the CMS study to assess the expected 
differences between the use of standard-cone isolation and a Frixione isolation for the pho- 



ton. To do so, we have made use of an NLO code, due to Gordon and Vogelsang [3l|. It is 
a semi-analytic code, relying on the narrow-cone approximation. We also used BlackHat 
in conjunction with SHERPA as a cross check on the NLO Frixione isolation result. (In all 
these calculations, oem was set to 1/137, for the reasons noted above.) 

The results of the comparison are shown in figs. [3] and [Has a function of the photon prp. 
Fig- El shows the CMS measurement, with its statistical and systematic uncertainties, along 
with the NLO predictions for a fixed-cone isolation as given using JetPhoX in ref. 1J|, and 
for the Frixione isolation criterion computed using the Gordon- Vogelsang code. The scale 
variation band shown was computed using JetPhoX. Fig. H] compares several different NLO 
predictions to the JetPhoX prediction as given in ref. ^14]: fixed-cone and Frixione isolation 
ones, using the Gordon- Vogelsang code, and one for the Frixione isolation, using Black- 
Hat with SHERPA. All the predictions here use the CTIO parton distributions j32|, which 
have as{Mz) = 0.1179. We see that the fixed-cone prediction from the Gordon- Vogelsang 
code, while not identical to the JetPhoX one, is in excellent agreement, well within the 
experimental systematic uncertainties. In addition, the predictions for the Frixione-style iso- 
lation are also quite close, within 2% at low pj, and within 1% at high p^. It is the large-p^ 
region that is relevant to the primary study we perform in this paper. We expect that the 
smallness of this difference will carry over to processes with multiple jets in addition to the 
photon present here. We will include this difference in our overall uncertainty estimate for 
the Z-to-photon ratios. (See ref. 33|] for a preliminary comparison between the two isolation 
criteria in the context of polarized pp scattering.) Fig. H] also verifies the agreement in results 
from the Gordon- Vogelsang code for the Frixione-style isolation with those from Black- 
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FIG. 3: A comparison of different NLO theoretical predictions for 7 + X production at the LHC 
at 7 TeV. The CMS data points p4] are shown (red) with combined experimental uncertainties; 
the prediction using using the Gordon- Vogelsang (GV) code and the Frixione isolation criterion is 
given by the dashed (blue) line; the JetPhoX prediction using a standard cone isolation is given 
by the solid (black) line. The lower panel shows the two theoretical predictions normalized to the 
CMS data along with the scale dependence (shaded gray) as determined using JetPhoX. 

Hat+SHERPA to within 1%, uniformly in In this particular BlackHat+SHERPA 
calculation, we impose no jet requirements, and impose a fixed maximum hadronic energy 
of 5 GeV in the isolation definition f l2.3p in place of Ej- e with e = 0.025 as done elsewhere in 
our 7 + 2-jet study. (The difference in cross sections between using a fixed energy of 5 GeV 
and the energy fraction used in our study is less than 1% for pj, > 50 GeV.) 

The rest of our study makes use of the MSTW08 parton distributions. Predictions with 
this set, which has a larger as{Mz) of 0.1201, are 4% higher than those shown in fig. [3l fairly 
uniformly in pj.; however, this difference is expected to become much smaller when taking 
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FIG. 4: A comparison of four different NLO predictions for 7 + X production at the LHC at 7 TeV. 
All predictions are normalized to the JetPhoX predictions shown in fig. [3l The prediction for a 
standard cone isolation using the Gordon- Vogelsang (GV) code is given by the dot-dashed (red) 
line; that for the Frixione isolation criterion using the GV code, by the dashed (blue) line; and 
that for the Frixione isolation using the BlackHat-|-Sherpa code, by the solid (black) line. The 
JetPhoX scale dependence band is shown in gray. 

the ratio of Z + 2- to 7 + 2-jet quantities. 



IV. PHOTON-TO-Z RATIOS AT THE LHC 



Our focus in this paper is on using distributions measured for 7 + 2-jet production to 
predict similar distributions assembled from missing Et + 2-jet events. Ideally, we would 
compute the ratio of fully-differential cross sections, 

da^/dA 



da-r/dA ' 

where 

da^ da^+^^''' 



dA dE]:' ^""'dri^'' j'^trf^i" ^""'dEr '"'dn 



2'"^ jet 7 J 2'"^ jet 



(4.1) 



(4.2) 



This information could be used experimentally by adjusting the weight of each 7 -|- 2 jet 
event, removing the photon, and then reweighting the sample to obtain unit- weight events. 
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This would provide a sample of estimated Z + 2 jet events with the Z decaying invisibly. 

With sufficient statistics we could effectively construct such a quantity using n-tuples, 
tracking the matrix-element weights and parton momenta. However, besides the practical 
issue of carrying out such a procedure at NLO, we would have no way to usefully present 
such a quantity in two-dimensional form. Accordingly, we will study the ratios of a variety of 
singly-differential distributions. This can be thought of as a projected version of the above 
ratio, with a projection onto one of the axes in the complete differential cross section. 



We use the anti-fcr jet algorithm [28| with clustering parameter R = 0.5, where R = 
{AyY + as usual in terms of rapidity and azimuthal angle. 

Our Monte Carlo set up allows us to study any infrared and coUinear-safe observable. 
In the present study, we focus on three sets of cuts, theoretical parallels to those used by 
the CMS collaboration [Ij and intended to be relevant in different regions of the SUSY 
parameter space. In order to display these cuts, we make use of a special definition of the 
total transverse energy, which we label H^^, as the sum of the transverse energies of all jets 
with pt > 50 GeV and |?7| < 2.5. We also define a vector MET, as the negative of the 
four-vector sum of all jets with px > 30 GeV and \ri\ < 5. Each set of cuts that we will 
consider is distinguished by different restrictions on the quantities^ Hj^^ and MET. 

Set 1: i/j.^* > 300 GeV, |MET| > 250 GeV, 
Set 2: if^"* > 500 GeV, |MET| > 150 GeV, 
Set 3: H^T^ > 300 GeV, |MET| > 150 GeV. 

For all sets we insist that the two highest-pr jets have pt at least 50 GeV and pseudora- 
pidity of at most |?7| = 2.5. These jets are referred to as 'tagging jets'. We note that beyond 
leading order there can be other jets in the event. The separation in 0-space between each 
tagging jet and the MET vector is required to satisfy A0(jet, MET) > 0.5. 

In addition to the above cuts, for the 7 -|- 2 jet study only, we impose photon isolation 



according to the Frixione 121] prescription, with parameters e = 0.025, 60 = 0.3 and n = 2, 



and a minimum i?-space separation between the MET vector and each tagging jet of 0.4. 



^ MET stands for 'Missing Transverse Energy'. Although this is standard terminology, it is potentially 
misleading, as the MET is in fact the missing transverse momentum. 
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The photon is required to have I77I < 2.5. We also impose a minimum pt of 100 GeV on both 
the photon and the Z. This cut has no effect because of the |MET| cut, but it improves the 
numerical efficiency of the calculation. 

The Set 1 cuts can be roughly characterized as the low-i/j?* / high-MET region, whereas 
Set 2 is the converse, high-i/j?* / low-MET. The reason for studying these two sets is that 
different SUSY production mechanisms are expected to lead to signals in different regions. 
Broadly speaking. Set 1 is geared towards catching direct squark decays, while Set 2 is 
designed for cascades with a 1^-boson and a softer lightest supersymmetric particle (LSP). 
Set 3, which is inclusive of both the others, is a control region. 

Our fixed-order results depend on the renormalization and factorization scales. These 
scales are unphysical, but necessarily appear when the pcrturbativc scries is truncated at 
a finite order. For fixed-order predictions, it is customary to estimate the error arising 
from omission of higher-order terms by varying these scales around some central value. 
The size of the resulting band is a useful diagnostic for those situations where fixed-order 
perturbation theory breaks down. The central value should be a typical hard scale in the 
process, to minimize the impact of potentially large logarithms. We choose the dynamical 
scale 11 = for this central value, where is defined as 

H!r = J2ET + ETiZ,7), (4.3) 

i 

where i runs over the partons and Et = \J M"^ + p\. A common method for estimating 
the error on a cross section is by varying the common scale up and down by a factor of 
two. We do so by evaluating the cross sections at five scales: /x/2, /x/-\/2, /x, -\/2a*, 2/x. As we 
will discuss below, this procedure is expected to greatly underestimate uncertainties when 
applied to a ratio of cross sections. 

V. LHC PREDICTIONS 

In this section we present total cross sections and distributions for 7-f 2-jet and Z+ 2-jet 
production at the LHC running at 7 TeV. We present results for each of the three sets of 
cuts. In the Z -\- 2 jet study we fold in the Z boson decay into neutrinos, which comprise the 
missing energy. The branching ratio for the Z to decay to neutrinos is largely responsible 
for the 7 -I- 2 jets cross section being about a factor of five larger than for Z{-^ vV) -\- 2 jets. 
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process 


LO 


ME+PS 


NLO 


Z + 2j 


0.521(0.001)1°;J^° 


0.416(0.004) 


0.560(0.002)l[!;j]^^ 


7 + 2i 


2.087(0.005)1|]:S^ 


1.943(0.027) 


2.448(0.008)l°:ii 


Z/7 ratio 


0.250 


0.214 


0.229 



TABLE I: Set 1 cross sections for Z and 7 production in association with two jets, using the anti-ZcT 
jet algorithm. The numbers in parentheses are Monte Carlo statistical errors, while the upper and 
lower limits represent scale dependence. See the text for a discussion of the errors on the ratio. 



process 


LO 


ME+PS 


NLO 


Z + 2j 


0.205(0.001)1^50 


0.238(0.005) 


0.277(0.002)l°:0i 


7 + 2i 


0.952(0.004)1[!;|3 


1.132(0.010) 


1.374(0.008)l[!:i| 


Z/7 ratio 


0.215 


0.211 


0.201 



TABLE II: As in Table 1, but for the Set 2 cuts. 



This ratio is clearly visible in our figures, and is of course the underlying motivation for this 
study. We will also discuss the error to be assigned to our predictions. 

In Tables [I| ITTl and UTTl we give the total cross section for the cuts outUned in the previous 
section. Each table shows three different theoretical predictions. Fixed-order perturbative 
results are shown as LO and NLO. The final states in these cases consist of the vector boson 
with the two tagging jets, though there can be an extra jet at NLO. The parton shower result, 
labeled "ME+PS" , is a tree-level matrix element calculation merged to a parton shower 11| . 
as summarized above. Here the final state can contain many jets, though virtual corrections 
are not taken into account. The LO predictions are the least reliable of the three and are 
shown only for reference purposes. 

With the cuts of Sets 1 and 2, the corrections from LO to NLO lead to an increase in 
the total cross sections of up to 50%. The corrections for the control region, given by the 
Set 3 cuts, are much more modest. The ME+PS and NLO results do not agree well for 
the cross sections. However, when one takes the ratio of Z and 7 cross sections, the two 
predictions agree to within better than 10%. This behavior is not surprising: typically, 
overall normalizations can be somewhat off in ME+PS calculations, while ratios tend to 
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process 


LO 


ME+PS 


NLO 


Z + 2j 


1.255(0.002)l|j:^^^ 


1.174(0.011) 


1.476(0.007)1[!;?^^ 


7 + 2i 


5.854(0.011)l}i6l 


6.075(0.056) 


7.601(0.019)10:81^ 


Z/7 ratio 


0.215 


0.194 


0.195 



TABLE III: As in Table 1, but for the Set 3 cuts. 



behave much better. 

We have found that in the ratio the LO scale variation cancels nearly completely, if we 
vary the scale in a correlated way in the Z + 2-jet and 7 + 2-jet predictions. In the NLO 
case the scale variation is a bit larger but also very small. This nearly complete cancellation 
of the scale variation cannot be interpreted as a small theoretical uncertainty. The closeness 
of the NLO and ME+PS ratios is a much better indication that the theoretical uncertainties 
on the individual cross sections do indeed largely cancel in the ratio. We do not include the 
uncertainty due to the parton distributions in our study, but we expect it to largely cancel 
in the ratio since the d{x)/u{x) ratio feeding into it is from a well-measured region in x. 

Let us turn next to an examination of five different distributions: the transverse momen- 
tum of the first (largest pr) accompanying jet; the pr of the second jet; the total transverse 
energy of the jets H^^^] and the azimuthal angle distributions with respect to the MET vec- 
tor of the first and second accompanying jets. These variables are useful for assessing the 
extent to which the kinematics of 7 + 2-jet events resembles that oi Z + 2-jet events. 

The pt distribution for the leading jet is shown in fig. |5l The peak in the distribution is 
well above the jet cut, because of the additional cuts on the MET and ifj^*. The LO cross 
section vanishes below 150 GeV for Set 1 and 250 GeV for Set 2, due to the restricted LO 
kinematics, which are relieved both at NLO, from the presence of one additional parton, 
or in the ME+PS approximation, which can have several additional partons (including the 
parton shower). Clearly fixed-order LO is inadequate for describing the pt distribution for 
the leading jet. 

The Pt distribution for the second accompanying jet is shown in fig. |6l we consider only 
Pt > 50 GeV. This distribution is peaked at the jet cut of 50 GeV for Set 1, and falls rapidly 
with increasing pt- For Set 2, in contrast, the distribution peaks at a second-jet pt of around 
180 GeV. The latter reflects a compensation for the higher total jet Et needed to arrive at 
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FIG. 5: The pT distribution of the first jet. The left column shows distributions for the Set 1 cuts, 
and the right column for the Set 2 cuts. Each column displays the differential cross section for 
Z + 2-jet production (top), 7+ 2-jet production (middle), and their ratio (bottom). In the top 
and middle plots, the upper panel shows the LO, NLO, and ME+PS results for the distribution, 
and the lower panel shows the ratio to the central NLO prediction, along with the LO and NLO 
scale-dependence bands. The numerical integration uncertainties are indicated by thin vertical 
lines. 1 7 
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FIG. 6: The distribution of the second jet. The plots are arranged and the curves are labeled 
as in fig. [5l 



the minimum ifj? in this set; a second jet at the lower cut would force the leading jet too 
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far out onto the tail of its distribution. The NLO corrections for Set 1 are approximately 
20%, in line with the corrections to the total cross section. These corrections are roughly 
flat across the distribution, so that the shape does not change at NLO. The NLO corrections 
for Set 2, in contrast, are significantly larger, up to 100% in some bins of second-jet pt, but 
under 10% for px above 300 GeV, but below 400 GeV, where the integration errors become 
large. The shape of the distribution correspondingly suffers significant corrections at NLO. 
For both sets, however, the corrections do not distinguish Z and 7 production, so that the 
ratio has NLO corrections of 10-15%, but fiat across all pxs. The NLO and ME+PS results 
are likewise in good and uniform agreement for the ratio, to better than 10%, again ignoring 
regions with large integration errors. As expected, the scale-dependence bands, from varying 
the common renormalization and factorization scale, = up and down by a factor of 
two, are narrower at NLO than at LO. 

The situation changes slightly if we examine the if^* distributions, shown in fig. [71 In 
both Sets 1 and 2, these are falling distributions which peak at the cut values for if;^* (300 
and 500 GeV, respectively). In Set 1, the NLO corrections range from 10% to 50%, gradually 
increasing with increasing H^^ through 550 GeV or so; the shape of the distribution is thus 
modified at NLO. In contrast, the corrections for Set 2 are fairly uniform across all H^x^, 
presumably because of the larger ifjT* starting scale for the distribution. Once again, the 
corrections do not distinguish between Z and 7 production, so that the ratios for both Sets 
1 and 2 have NLO corrections of about 10-15%, with corrections relatively fiat across all 
H'^\ The NLO and ME+PS results for the ratio are again in good and uniform agreement, 
to better than 10%. In both sets, the scale-dependence bands are uniformly narrower at 
NLO than at LO. 

The distributions which show the most structure are those in the azimuthal angle 
A0(MET, jetj) between the MET vector and the two jets, shown in fig. [HI The distribution 
for the first jet is peaked near A0 = tt, where the leading jet balances the MET vector while 
the second jet has much smaller pt- The distribution falls very steeply at smaller angle. A 



typical LO configuration for the Set 1 cuts is shown in fig. 9(a) With only two partons and a 
missing energy vector from the Z or 7, transverse momentum conservation implies that the 
leading jet can never get closer in angle than 7r/2 to the MET vector. It can only approach 
this angle if the transverse momenta of both jets are extremely large. This configuration is 
heavily suppressed by the parton distributions, to orders of magnitude below what is shown 
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FIG. 7: The Hj^ distribution. The left column shows distributions for the Set 1 cuts, and the 
right column for the Set 2 cuts. The plots are arranged and the curves are labeled as in figs. [5] and 
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FIG. 8: The Ai;^(MET, jet) distributions for the leading two jets. The left column shows the 
distributions for the Set 1 cuts, and the right column for the Set 2 cuts. Each column displays 
the differential cross section for Z + 2-jet production (top), 7+ 2-jet production (middle), and 
their ratio (bottom). The left panel of each plot shows the distribution for the leading jet, and the 
right panel for the second jet. In the top and middle plots, the upper panel shows the LO, NLO, 
and ME+PS results for the distribution, and the lower panel shows the ratio to the central NLO 
prediction, along with the LO and NLO scale-dependence bands. 
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FIG. 9: A typical configuration of jets and MET vector for the Set 1 cuts at (a) LO and (b) NLO. 

in fig. [HI At NLO, this kinematic constraint is relaxed, as it is possible in the real-emission 
configurations for the second and third jets (the latter may fall below the jet px threshold) to 



balance the leading jet and MET vector, as shown schematically in fig. 9(b), At this order, 
the leading-jet distribution does indeed run all the way to the A0(MET,jet) = 0.5 cut we 
impose. The cross section for angles less than 7r/2 is tiny, however, so although the NLO 
corrections are very large in this region, they are of no practical importance; this carries 
over to the large corrections seen in the ratio of Z + 2-jet and 7 + 2-jet distributions. 

The width of the peak region in the A0 distribution for the leading jet can be understood 
by considering a configuration in which the second jet has pr — 50 GeV and 180 GeV for 
Sets 1 and 2 respectively, with the leading jet having the minimum p^ needed to balance 
the MET given a jet-jet A0 = 0.5 (the interjet AR cut). This would lead to a peak around 
A(f) 3.05 for Set 1, and 3.0 for Set 2. However, the shapes of the ratios of distributions do 
differ somewhat between the NLO and ME+PS predictions for the peak region of the Set 1 
distribution (A0 ~ 2.8-3.0), so that the analysis would be more sensitive to any systematic 
effects in the treatment of this region. 

In the A0 distribution for the second jet, we see different structure in the small-angle 
region for Z + 2-jet and 7 + 2-jet production. For Set 1, this region is subdominant in 
its contributions, although not completely negligible, because the high MET cut favors 
configurations in which both jet transverse momenta oppose the MET vector, as shown 
in fig. [Ht^a). In Set 2, however, the low MET cut allows small values of A0(MET,jet2) 
to dominate. In this region, the potential collinear singularity between a photon and an 
outgoing quark forming the second jet causes the cross section to start to rise at small 
angles. This rise is cut off by the A0(MET,jet) cut, but does distinguish Z production 
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from photon production. For both Sets 1 and 2, the shapes of the ratios of distributions 
are similar at LO, at NLO, and in the ME+PS resuhs, so that we can be confident the 
Z-to-7 difference is captured properly, and that the analysis should not be overly sensitive 
to the precise value chosen for the A0(MET, jet) cut. The shapes of the NLO and ME+PS 
distributions (though not the LO one) are quite similar as well in the peak region, so that 
the analysis should be robust under systematic effects treating different angular regions in 
A0(MET,jet2) somewhat differently. 

To assess the overall uncertainty, we see that the ME+PS prediction for the total cross 
section with the Set 1 cuts is approximately 7% lower than the NLO prediction, while for 
the Set 2 cuts it is approximately 5% higher. In addition, we should allow for a difference 
of up to 1% in the prediction for the fixed-cone isolation compared to the Frixione- style 
isolation, as discussed in section lllli This suggests that the NLO predictions for the ratios 
used by the CMS collaboration should be accurate to within 10%. 

We have not considered the effect of possible electroweak Sudakov logarithms. Based on 
refs. {341, we estimate that for Set 1, the effect will be under 5% and smaller for Sets 2 and 
3. Our study considers production of an (exclusive) electroweak boson accompanied by at 
least two jets rather than one; we expect the electroweak effects to be of similar magnitude 
for similar boson transverse momenta. If one were to raise the MET cut, these electroweak 
corrections increase, and would need to be taken into account to ensure the reliability of 
the prediction for the Z + 2-jet to 7 + 2-jet ratio. Were the MET cut to be raised to 
500 GeV, for example, the virtual electroweak corrections would grow to about 10%. The 
effects of radiating a second electroweak gauge boson into the final state |35i] should also be 
considered. These questions merit further investigation. 



VI. CONCLUSIONS AND OUTLOOK 



Data-driven estimates of backgrounds to new physics signals offer a powerful means for 
avoiding a reliance on Monte Carlo estimates of backgrounds. Even in such data-driven 
methods, theoretical input is usually required to provide conversion factors from one process 
to another, and estimates of uncertainties in these factors, as well as of extrapolations from 
control to signal regions. 

In this paper, we studied the theoretical aspects of using 7+ jets data to estimate the 
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missing- i?T + jets background to new-physics searches. We focused on the comparison of 
boson production in association with two or more jets, computing the relevant differential 
cross sections and ratios in NLO QCD. Because the scale uncertainty is quite small in these 
ratios, we used the difference between ME+PS results and the NLO ones as an estimate 
of theoretical uncertainty. (We did not study the uncertainty due to parton distribution 
functions, but it is expected to be small as well.) Our study used the Frixione isolation 
criterion to compute the prompt-photon cross sections; a comparison with the fixed-cone 
isolation for isolated prompt-photon production indicates that the resulting shift should be 
less than 1% in the high-pj. region of interest. We did not include the effects of electroweak 
Sudakov logarithms, and while these will become significant for higher MET cuts, we expect 
the effects to remain small for the cuts used in this study. 

We find that the conversion between photons and Z bosons has less than a 10% theoretical 
uncertainty for events with two associated jets. In the future, it should be feasible to extend 
this study to three associated jets. The small uncertainty we find should make it possible for 
the photon channel to provide a competitive determination of the Standard-Model missing- 
Et + jets background for many years to come. 

Note added 

Since the appearance of our paper, CMS presented a new physics search based on an event 



signature of at least three jets accompanied by large missing transverse momentum [36|. No 
excess events were observed above the background. The irreducible background from Z 
bosons decaying into a vv pair was estimated by converting a measurement of photons 
accompanied by jets into a prediction of this background. Our paper provided the estimate 
of the theoretical uncertainty for this conversion. 
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